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SUMM^Y 

In  this  study  a mathematical  method  to  determine  the  response  of  the  blood  vessels  to  transient  and 
sustained  acceleration  forces  is  presented.  The  method  is  based  on  coupling  of  the  Navier-Stokes  equations 
for  blood  flow  and  the  large  elastic  deformation  theory  for  the  deformation  of  the  blood  vessels,  and 
solving  them  numerically  under  the  appropriate  initial  and  boundary  conditions.  A mathematical  reasoning 
to  neglect  the  effect  of  acceleration  on  microcirculation  per  ae  is  given.  However,  microcirculation  is 
indirectly  affected  by  acceleration  forces  which  tend  to  pool  blood  and  bring  about  pressure  changes  in 
large  vessels.  Aortic  pressures  are  calculated  for  examples  of  monotonically  increasing  and  transient  -G 
acceleration  profiles,  and  one  of  the  solutions  is  compared  with  an  available,  experimentally  measured 
pressure  from  an  animal  experiment.  In  the  absence  of  proper  physiological  scaling  laws,  the  qualitative 
agreement  between  the  theory  and  experiment  is  satisfactory.  > 


INTRODUCTION 

Recent  developments  in  spacecraft  and  high  performance  aircraft  have  resulted  in  the  exposure  of  the 
human  body  to  the  hazards  of  high  accelerations  beyond  tolerance  levela.  In  addition  to  injury  to  various 
parts  of  the  body,  cardiac  insufficiency  and  the  consequent  physiological  malfunctions  such  as  headache, 
abdominal  pain,  change  in  heart  rate,  chest  pain,  loss  of  vision  and  hemorrhage  are  seme  of  the  manifesta- 
tions of  acceleration  trausm.  The  knowledge  of  the  response  of  the  cardiovascular  system  to  acceleration 
stress  is  essential  to  the  development  of  protective  devices  which  are  designed  to  increase  acceleration 
tolerance  by  the  human  body  during  aircraft  and  spacecraft  maneuvers.  The  objective  of  the  present  in- 
vestigation is  to  understand  the  response  of  the  blood  vessels  to  transient  and  sustained  aerospace  accel- 
eration profiles. 

Since  it  is  impossible  to  actually  subject  the  human  body  to  abnormal  high  accelerations  without  in- 
flicting an  injury,  the  response  of  the  vascular  system  to  transient  and  sustained  accelerations  is  investi- 
gated mathaewtically.  Theoretical  analyses  are  extremely  helpful  for  evaluating  the  relative  injury 
potential  for  various  acceleration  functions,  in  guiding  experimental  investigations,  and  in  developing 
and  understanding  protective  measures.  Mathmaatical  procedures  also  provide  the  basis  for  establishing 
precise  dynamic  and  physiological  scaling  laws  needed  to  translate  experimental  data  obtained  with  various 
species  into  meaningful  results  for  himans. 

In  recent  years,  numerous  mathematical  investigations  of  arterial  dynamics  have  appeared  in  the 
scientific  literature.  Momeraley  11]  and  Noordergraaf  U)  have  presented  a mathematical  analysis  of  blood 
flow  through  arteries  by  using  a lumped  parameter  model.  Taylor  13),  Kenner  14),  Attinger  et.  al.  (5)  used 
distributed  parameter  models  to  analyze  pressure-flow  relationships  in  arteries  and  veins.  Several  articles 
related  to  blood  flow  in  arteries  have  appeared  in  the  book  by  McDonald  [6) . An  elastic  tube  theory  of 
blood  flow  has  been  treated  by  Lambert  17]  and  Skalak  and  Stathis  [8] . Kivity  and  Collins  19]  presented 
a viscoelastic  tube  model  for  aortic  rupture  under  decelerative  forces.  Rudinger  110)  studied  the  effect 
of  shock  waves  on  mathematical  models  of  aorta  for  better  understanding  of  the  behavior  of  the  actual  aorta. 

To  understand  the  blood  flow  characteristics  in  the  arterial  system,  the  knowledge  of  the  material 
properties  of  the  arterial  wall  is  essential.  Bergel  111],  Fung  [12],  Demiray  and  Vito  [13]  have  utilized 
matheamtical  models  of  the  constitutive  proparties  of  the  arterial  tissue  to  determine  the  stresses  in  the 
arterial  walls.  In  the  present  study,  the  strain  energy  function  given  by  Demiray  and  Vito  [13]  for  an 
arterial  wall  specimen  has  been  used  in  determining  the  aortic  pressure  that  is  compatible  with  large  de- 
formation of  the  aorta  and  the  associated  flow  under  acceleration  stress. 

The  present  study  also  considers  the  effect  of  acceleration  on  the  microcirculation.  Microcirculation 
under  normal  conditions  wee  investigated  by  Pro the ro  and  Burton  (14),  Whitmore  [IS],  Gross  and  Aroesty  (16), 
Gross  and  Intaglietta  [17] , Skalak  [18]  and  Fung  [19]  who  presented  various  theories  of  flow  in  the  capil- 
lary bed  connecting  the  arteries  and  veins. 

Several  axperlawntal  investigations  on  the  affects  of  acceleration  stress  on  the  human  body  have  bean 
performed  at  the  USAF  School  of  Aerospace  Medicine  at  Brooke  Air  Force  Baas,  Texas,  fcirton  [20]  eubjzcted 
miniature  swine  to  Gs  acceleration  to  study  its  effects  on  the  organism  and  extrapolated  the  results  to 
human  beings.  Parkhurst,  et.  al.  [21]  conducted  experiments  on  human  tolerance  to  high  +C,  forces. 

Leverett,  et.  al.  [22]  investigated  the  physiologic  response  to  high  mistaimad  acceleration  stress.  Peter- 
son, et.  al.  [23]  studied  the  cardiovascular  responses  during  and  following  exposure  to  +G,  foroes  in 
chronically  lnatnxNnted  anesthetized  doge.  Burton  and  MacKensle  (24)  determined  the  extent  of  heart 
pathology  as  a function  of  acceleration  stress. 
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MATHEMATICAL  FORMULATION 

A.  Equations  of  Fluid  Motion 

The  geometry  of  the  elastic  tube  containing  blood  In  notion  is  shorn  in  Fig.  1.  Let  r,  6,  z be 
the  cylindrical  polar  coordinates  and  let  u,  v,  and  w be  the  velocity  components  in  the  corresponding 
directions.  Assuming  axial  synaetry  in  flow  and  tube  deformation,  the  Navier-Stokes  equations  for  the 
flow  of  blood  can  be  written  as : 


3u  3u  3u  1 3p  . r a2u  . 1 3u  32u  u 1 

37+u37  + w37"-7;  a?+vL»?2  r»?  + II?-7rJ 

3w  3w  3w  1 to  f 3^  1 3w  3^  "1  ... 

3?  + u Jr  + w 7s  ” ~ 3s  + v * 7 3r  + 3*1"  J * 9(t) 


where  p is  the  pressure,  v is  the  kinematic  viscosity,  pQ  is  density  of  blood  and  g(t)  is  the  body  force 
per  unit  mass  caused  by  the  acceleration.  The  continuity  equation  is 

2“  ♦»+!«.  0 (3) 

3r  r 3z 

The  above  equations  are  nondimensionalised  using  a typical  length,  Rq,  which  is  the  initial  (undeformed) 
radius  of  the  aorta,  and  U,  the  average  velocity  of  blood  in  the  aorta.  Introducing  the  new  quantities 
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the  equations  of  motion  and  the  continuity  equation  in  terms  of  the  newly  defined  variables  become 
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Deleting  the  "stars”  for  simplicity,  the  governing  equations  in  the  dimensionless  form  will  becosm 
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The  boundary  and  initial  conditions  are 
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t » 0 

where  is  the  inside  radius  of  the  blood  vessel  in  the  deformed  state. 

B.  Equations  of  Notion  for  Thin-welled  Elastic  Tubei 

The  theory  of  large  elastic  deformations  is  utilised  to  describe  the  time-dependent  deformation 
of  the  blood  vessels . In  view  of  the  published  results  on  blood  pooling  and  the  consequent  cardiac  in- 
sufficiency, the  application  of  large  deformation  theory  appears  necessary.  Damiray  and  Vito  [13]  have 
previously  used  this  theory  to  calculate  the  deformation  of  arteries. 

The  undaformed  and  deformed  cylindrical  tubes  are  shown  in  Fig.  2.  Let  r,6,  s represent  a point  in 
the  wall  of  the  undeformed  tube,  and  X,  6,  i in  the  deformed  tube,  r,,  r are  inside  and  outside  radii, 
respectively,  of  the  undeformed  tuba,  and  R. , R those  of  the  deformed  tuba.  Axial  stretch  of  the  tube  is 
neglected  because  of  tethering  caused  by  the  surrounding  tissue.  As  evening  the  material  of  the  blood  vessels 
to  be  homogeneous,  incompressible,  and  isotropic,  the  stress  at  any  point  can  be  written  aai 

■ eg11  ♦ *BlJ  ♦ PC11  <M> 


4 ■ 2(9W/3l^) , 4 ■ 2(3W/?I2)  , B J ■ I^g  J - 9 gJ  Gra#  P is  a scalar  function  which  represents  a 
hydrostatic  pressure#  N is  the  strain  energy  function#  I±  and  12  are  the  strain  invariants#  and  g*J#  g^j# 
G*3,  and  G^j  are  the  contr avar iant  and  covariant  netric  tensors  (25,  26].  The  indices  i and  j take  the 
values  1#  2#  and  3.  The  equations  of  notion  are  given  by: 

riJll  i + PwrJ  - P„f3  (13) 

where  |)  denotes  covariant  differentiation,  p„  la  the  density  of  the  vessel  vail,  P Is  the  body  force,  and 
f Is  the  acceleration.  Let  us  neglect  the  body  force  on  the  vessel  wall  In  comparison  to  its  effect  on 
the  fluid  flowing  in  the  cylindrical  tube.  Performing  the  covariant  differentiation  on  tne  remaining  part 
of  the  equation  of  motion  we  get 


in  which  a and  6 are  material  constants.  Defining  the  circumferential  stretch  ratio  X 
in  the  r,  8,  s directions  can  be  expressed  as 
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The  initial  conditions  are: 

At  time  t - t^#  R ,•  R , dR./dt  - u,  radial  velocity  of  fluid, 
o 1 o 1 


C.  Effect  of  Acceleration  on  Microcirculation 

The  blood  vessels  of  microcirculation  are  extraordinarily  small*  and  their  typical  dimensions 
are  of  the  order  of  microns.  Under  normal  circumstances#  the  velocity  of  the  blood  in  the  microcirculation 
is  1 mm/ sec  and  the  Reynolds  number  is  of  the  order  0<1<T3)  # which  is  sufficiently  small  so  that  the  Stokes 
flow  approximations  are  applicable.  Neglecting  the  inertial  effects  and  assuming  that  the  stream  lines 
are  nearly  parallel#  the  dimensionless  equation  of  fluid  motion  in  the  axial  (z)  direction  becomes 
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which  can  be  rearranged  to  read 

-If — (15  “ 12,1 

In  the  earth's  natural  gravitational  field,  the  dimensionless  g.  as  given  in  Eq.  (4).  is  of  the  order 
0(10"2),  and  with  the  effect  of  Re  - O(10*3)  in  the  last  tern  Re  g(t)  in  Eq.  (27)  becones  physiologically 
insignificant,  being  of  the  order  0<10_s) . Re  estimate  that  the  effect  of  acceleration  on  microcirculation 
per  se  can  be  safely  neglected  up  to  100  g.  However,  the  pressure  of  the  blood  pooled  in  the  arteries  and 
veins  can  affect  the  flow  rate  in  the  nail  vessels.  For  this  reason  it  is  necessary  to  determine  a re- 
lationship between  the  pressure  gradient  and  the  flow  rate  in  the  small  blood  vessels. 

For  the  flow  of  a Newtonian  fluid  in  a uniform  tube  Szymanski  (27)  showed  that  the  flow  would  be 
fully  developed  if  vt/D2  V 1,  where  t - time,  v - kinematic  viscosity,  and  D - tube  diameter.  An  extension 
of  this  criterion  to  microcirculation  yields  vAt/D2  V 1 for  flow  to  be  quasi-steady,  where  At  is  the 
nallest  characteristic  time  of  the  unsteadiness  in  flow.  According  to  Burton  (28),  At  or  0.1  sec;  using 
v « 0.04  Stokes,  one  finds  that  the  diameter  D must  be  greater  than  600(j  (microns)  for  any  significant 
effect  of  unsteadiness.  Since,  in  microcirculation  the  diameters  of  blood  vessels  are  much  less  than  600u, 
changes  in  flow  due  to  unsteadiness  become  entirely  negligible.  On  this  basis  Benis  (29)  argued  that  the 
effect  of  unsteadiness  on  non-Newtonian  flow  could  also  be  neglected.  Thus,  the  use  of  steady-flow  equa- 
tions can  be  justified  for  microcirculation. 

For  steady  capillary  flow,  the  flow  rate  through  a circular  tube  can  be  expressed  by 


0 - 2w  j rw  dr  (28) 

o 

where  Q - flowrate,  R - tube  radius,  and  w - blood  velocity.  Integration  by  parts  of  the  right  hand  side 
yields 


**(£)-'  ,29) 

The  first  integral  on  the  right  hand  side  of  Eq.  (29)  is  sero.  In  the  second  integral  the  domain  of  in- 
tegration can  be  divided  into  two  regions!  a cone  of  unshearad  fluid  extending  to  radius  Ry,  and  the 
annular  region  bounded  by  the  unsheared  fluid  and  the  tube  wall. 

Than. 

*y 

The  first  term  on  the  right  hand  side  in  Eq.  (JO)  which  represents  the  core  integral  is  sero.  By  changing 
the  variables  from  r to  l in  the  second  term  as  suggested  by  Merrill  et.  el.  (30),  Eq.  (30)  can  be  written 


»-  T&ST  /■’■*  * 

T 

y 

where  L - length  of  the  capillary 
AP  ■ pressure  drop 
t • shear  stress 
f » sheer  rate. 


The  shear  stress  and  the  shear  rate  are  related  by  an  empirical  equation 


T 1/2  . tyl/2  + yl/2  y 1/2 


(32) 


In  which  Ty  Is  the  yield  shear  stress,  and  y'^2  is  a constant  which  represents  the  slope  of  the  Casson 
plot  relating  the  viscometric  parameters  of  blood.  In  Foiseuille's  flow  y becomes  the  blood  viscosity. 
Substitution  of  Eq.  (32)  Into  Eq.  (31)  and  Integration  yields 
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which  is  valid  under  the  assumption  that  the  flow  is  steady,  laminar  and  incompressible,  and  blood  is 
homogeneous.  In  the  above  equation,  Ty  and  y are  known  constants)  then  plots  of  Q vs.  AP/L  for  capillaries 
of  various  radii  can  be  easily  constructed.  An  example  of  this  relationship  is  shown  in  Fig.  6. 


RESULTS  AND  DISCUSSION 

To  determine  the  response  of  the  blood  vessels  one  must  solve  Eqs.  (8),  (9),  (10)  and  (25)  under  the 
appropriate  boundary  conditions.  Since  these  equations  were  coupled  a numerical  solution  using  a digital 
computer  was  sought  as  described  in  an  earlier  paper  [31] . 

The  following  constants  were  used  in  the  solution:  1^  - 1.47  cm,  0 - 11.9  cm/s,  pw  - 1.05  gr/cm3, 

PQ  • 1.05  gr/cm3,  v - 0.038  Stoke,  a * 0.8  and  8 - 11.35  x 104  dynes/cm2.  The  constants  a and  6 are 
elastic  constants  which  appear  in  the  strain  energy  function  W,  Eq.  (15) . These  values  are  reported  in 
literature  for  a specimen  of  htmum  aorta.  The  elastic  constants  for  specimens  of  veins  are  not  available, 
primarily  because  the  research  of  determining  the  strain  energy  function  in  the  form  of  Eq.  (15)  is  rela- 
tively new. 

As  examples,  two  deceleration  profiles  were  used  in  the  solution:  one,  a linear,  monotonieally  in- 
creasing function  of  time  represented  by  g(t)  ■ 7840t  cm/s2,  and  a transient  type  (Fig.  3)  which  increased 
and  decreased  rapidly  [32].  For  both  these  cases,  aortic  pressures  were  calculated  by  utilising  a finite 
difference  technique  which  involved  Runge-Kutta  integration  procedure  and  Adau-Bashforth  predictor- 
corrector  method.  The  dimensions  of  the  human  aorta  were  chosen  frem  the  physiological  data  presented  in 
Westerhof  et.  al.  [33].  The  aortic  pressures  calculated  for  the  sustaining  linear  deceleration  profile 
is  shown  in  Fig.  4.  The  aortic  pressure  in  response  to  the  transient  profile  of  Fig.  3 is  presented  in 
Fig.  5 along  with  an  experimentally  determined  pressure  in  the  thoracic  aorta  of  a beagle  dog  which  was 
subjected  to  the  same  transient  deceleration  profile  for  comparison.  The  shapes  of  the  pressure  vs.  time 
curves  are  nearly  the  same  indicating  satisfactory  qualitative  agreement.  An  exact  quantitative  agreement 
cannot  be  expected  because  of  the  anatomical  and  physical  differences  between  the  subjects. 

The  method  of  calculation  described  above  for  an  artery  can  be  easily  extended  to  a vein,  and  the 
pressure  difference  across  the  capillary  bed  can  be  determined.  Equation  (33)  represents  a relationship 
between  the  flow  rate  and  the  pressure  gradient  in  a small  tube.  For  selected  values  of  ty  “ 0.042 
dyne/cm2,  y “ 0.05  dyne-sec/cm2,  D - 100  microns,  L ” 2.5  cm,  the  flow  rate  Q is  computed  for  various 
pressure  drops  in  the  range  AP/L  » 10  —10000  dyne/cm3  and  shown  in  Fig.  6.  If  the  pressure  of  the  blood 
pooled  on  the  venous  side  is  known,  the  effect  of  acceleration  on  the  flow  in  microcirculation  can  be 
determined  from  Fig.  6.  It  must  be  noted  that  the  flow  rate  indicated  in  Fig.  6 is  for  a narrow  blood 
vessel  of  specific  dimensions.  To  obtain  the  total  blood  flow  one  must  formulate  the  solution  on  a 
statistical  basis  which  includes  the  arterioles,  venules  and  capillaries  of  various  dimensions  and  chang- 
ing rheological  properties  of  blood. 

The  results  presented  in  this  paper  are  part  of  an  effort  to  describe  the  cardiovascular  system  in 
terms  of  its  physical  and  mechanical  properties.  Most  investigations  hither  to  reported  have  dealt  with 
electrical  analogs  of  the  cardiovascular  system  in  which  various  parameters  were  introduced  in  terms  of 
resistances,  impedances  and  capacitances.  These  electrical  quantities  may  not  truly  represent  the  cardio- 
vascular parameters  under  high  g conditions,  and  an  analysis  of  the  cardiovascular  systM  which  is  based 
on  its  original  properties  is  therefore  desirable. 
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